COMMUTATIVE PARTIAL DIFFERENTIAL OPERATORS 



ALEX KASMAN AND EMMA PREVIATO 

Abstract. In one variable, there exists a satisfactory classification of commutative 
rings of differential operators. In several variables, even the simplest generalizations 
seem to be unknown and in this report we give examples and pose questions that 
may suggest a theory to be developed. In particular, we address the existence of 
a "spectral variety" generalizing the spectral curve of the one dimensional theory 
and the role of the differential resultant. 



1. Introduction 

In one variable, there exists a satisfactory classification of rings of differential op- 
erators that are (maximal) commutative. In several variables, even the simplest 
generalizations seem to be unknown and in this report we give examples and pose 
questions that may suggest a theory to be developed. To motivate what we do, we 
briefly recall the 1-variable case and selected results in several variables. 

In the 1-variables case with analytic coefficients the classification was found by 
Burchnall and Chaundy in the 1920s ]3J, [5j by essentially formal methods of differential 
algebra and some algebraic function theory on the "spectral curve" Spec A, where A 
is the commutative ring, completed by a point "at infinity". The isospectral rings, 
roughlty speaking, form the Jacobi variety of this projective algebraic curve and 



the "Krichever map" |1J| solves the inverse spectral problem explicitly. For a ring 
A generated by a pair L = d n + w n _2<9 n_2 + • • ■ + Uq{x) (where d = d/dx) and 
B = d m + ■ • • the determinant of the (n + m) x {n + m) resultant matrix of L — A and 
B — fi [2]J is a non-zero, polynomial p £ C[A, n] such that p(L, B) = 0. Moreover, in 



the "rank 1 case" gcd(n, m) = 1 |2l] one has that p(A,/i) = is the affine equation 
of the spectral curve. A divisor of points on the curve corresponding to the given 
element of the Jacobi variety is the set of poles of the (normalized) gcd of L — A and 
B — fi at a reference point x = xo; the flow in x is linear on the Jacobian. 

In view of this, it is natural to aks at least the following questions in several 
variables: 

1. Is Spec A an affine variety of dimension N for any maximal commutative ring 
A of partial differential operators in N variables? 

2. For a ring A with generators L,i (1 < i < N + 1), what is the relationship 
between the differential resultant of Li — Aj and the equation of such a variety? 
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More ambitiously, of course, one would ask for isospectral flows, a compete classi- 
fication, and the inverse spectral problem, in increasing order of magnitude. A beau- 



tiful generalization of Burchnall-Chaundy theory was given by Nakayashiki [18, 19 
using the Fourier-Mukai transform. He associates commutative rings in N variables 
to a suitable iV-dimensional abelian variety an some additional choices, one for each 
element of its Picard variety. But these are not scalar operators, rather they have 
(iV! x N\) matrix coefficients. 

In this report, we give a negative answer to question (1) by using techniques de- 
veloped in jnj] and we offer some observations, natural conjectures, and a strategy 
to treat (2). We deal with the scalar case only. We include as an appendix the Math- 
ematica code that can be used to compute the differential resultant of any set of 
N + 1 partial differential operators in N variables, which is a handy tool for checking 
properties on the available class of examples such as 0,0. 



2. Geometric Structure of Maximal Commutative Rings 

It is well known that the commutative rings of ordinary differential operators are 
finitely generated rings coordinate rings for algebraic curves. This forms the foun- 
dation of the Burchnall-Chaundy theory of such rings |4], |J |2l|. In contrast, very 
little is known about the algebro-geometric structure of commutative rings of partial 
differential operators. A question of interest is to address the problem of whether 
every commutative subring is contained in a (larger) commutative ring requiring only 
a finite number of generators over C. Such a result is relevant, for instance, to the 
algebro-geometric investigations of quantum integrable systems || . 

In this section we use techniques from soliton theory (namely Darboux transfor- 
mation and Baker- Akhiezer functions) to study the structure of certain commutative 
rings of differential operators. We are able to show that these rings are maximal in 
the sense that they are not contained in any larger commutative subrings of the ring 
of differential operators. This in itself is a difficult task which is rarely achieved. 
Then, in one example, we study the structure of this ring more closely and note that 
it cannot be constructed with only a finite number of generators over C and hence is 
not the coordinate ring of an affine algebraic variety. 

Notation: Let T> = C(xi, • • • ,x n )[d\,... ,d n ] be the ring of rational coefficient 
differential operators in n variables. It will be useful to be able to refer also to 
V = C[di, . . . ,d n ] C V (the constant coefficient differential operators) as well as 
a ring of microdifferential operators containing T> as well as the inverse of the 



particular operator K e T> which will be important below [JT3J and ^T> , the contant 
coefficient microdifferential operators. 



2.1. What Commutes with Many Constant Coefficient Operators? In gen- 
eral, it is difficult to address the question of whether a commutative ring of par- 
tial differential operators is maximal. The key which allows us to do it here is the 
observation that although the centralizer of a single constant coefficient operator 
p(d\, . . . ,d n ) will contain non-constant coefficient operators, an operator commutes 
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all constant coefficient multiples of p(di, . . . ,d n ) if and only if it also has constant 
coefficients. 

Lemma 2.1. Let p G TJ be a non-zero constant coefficient differential operator. 
Then any operator L G \I/T> which commutes with p as well as all operators Pi : = 
di ap(d\, . . . ,d n ) (1 < i < n) is also constant coefficient (i.e. L G ^T> ). 

Proof. Let us suppose that [p, L] = [Pi, L] = 0. Then since p o L = Lop and 
p o di = di o p it follows that 

= [P^ L] = di o p o L — L o diO p 
= di o L o p — L opo di 
= [Lop, di] 

Then we note that letting Lop = fa{x\, . . . , x^di 1 ■ ■ ■ d%™ be a series represen- 
tation for any microdifferential operator then 

[d l ,L\ = Y J a^... , Xn )d? i ---d? 

where prime denotes differentiation with respect to Xj. Hence, if L o p commutes 
with each di then Lop has constant coefficients. However, this provides a linear 
relation between any coefficient of L and certain higher ones. Since the coefficients 
of L are bounded factorially with the order, this is only possible if the coefficients are 
all constant. □ 



2.2. A Maximal Commutative Ring from Darboux Transformation. Here 
we will consider a special class of commutative rings of differential operators for which 
we are able to demonstrate maximality using the results of the previous subsection. 
Suppose that the constant coefficient operator p ET> factors as 

(1) p(d l} ... ,d n ) = LoK L,K <EV. 



(See and \T%[ for a discussion of some methods for achieving such factorizations.) 
Then the method of Darboux transformation commonly used in the study of inte- 
grable systems [Jj, ^ [| [H]] is to consider the (more complicated) operator 

P ■= K oL = K opoR- 1 

which shares many features with p since the two operators are conjugate. For in- 
stance, one may try to conjugate other constant coefficient operators by K to pro- 
duce operators that commute with P. In fact, given any constant coefficient oper- 
ator r(di, . . . ,d n ) G V it follows that [K o r o K~ l ,P\ = . However, although 
K o r o K^ 1 G ^T> there is no reason to expect that it is in T>. The content of the 
next theorem is the statment that the ring of all differential operators commuting 
with P which are produced in this way is a maximal commutative ring. 

Notation: Let K G T> be a differential operator and define 

R (K) :={reV \Ko r(d h . . . , d n ) G V o K} C P 
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to be the subring of elements r £ T>q such that K or has K as a right factor. Then 
the ring 

:= (iToDoo^jnD 

= k o Rq(k) o fsr 1 

is a commutative subring of P. In general, it will be the case that R(K) = C is 
trivial, but if K is chosen to be a non-constant operator satisfying (P then R(K) 
will contain differential operators. 

Theorem 2.1. Let K £ D be a differential operator which is the right factor of some 
constant coefficient operator p = L o K £ V . Then if R' is a commutative ring such 
that 

R{K) C R! C V 

it follows that R{K) = R' . In other words, R{K) is a maximal commutative subring 
ofV. 

Proof. Let us suppose that Q £ T> commutes with every element of R(K). We must 
show that Q is already in R(K). Note that since p = L o K one automatically has 
that Pi := di op is in Rq{K). Thus, we know that 

[K o P { o K~\ Q] = for 1 < i < n 

and hence conjugation in the ring of pseudo-differential operators gives 

\PuK~ 1 oQoK] = l<i<n. 



By Lemma 2.1, this implies that 



R- 1 o Q o K £ ^X?o 

and so it is clear that Q is of the form K o q(di, . . . , d n ) o K^ 1 for some constant 
coefficient pseudo-differential operator q £ ^T> . 

However, we may moreover note that q £ T>$ C ^T>o is a constant coefficient 
differential operator. To show this we introduce a normalized common eigenfunction 

... ,x n ,z u ... ,z n ) = -K [ e «i*i+-+*n*»] 

gr^l, . . . , Z n ) 

where the polynomial g is to be defined below. Note that regardless of the choice of 
g, we have by construction that 

K o L[ip) = p(z 1 , ... , z n )ip Q[ip) = q(z lt ... , z n )ip. 

One can write the function i^[exp(^ XiZi)] in the form 



j£ |" e xi«i+-+a5n«n] _ I ^a=l Pa( Z l' ■ ■ ■ I ^n)c r a(^l, ■ ■ ■ , X n ) i ^ 

' eroOi, ... ,x n ) ' 



where p a are all non-zero polynomials in zi, . . . ,z n and a a are distinct, non-zero 
monomials in xi, . . . , x n . We choose . . . ,z n ) to be the highest common fac- 
tor of the polynomials p a (zi, ■ ■ ■ ,z n ). We have thus constructed ip so the prod- 
uct f(zi, . . . ,z n )ip is holomorphic in each Zi for / £ C((z±, ... ,z n ) if and only if 
/ £ C[zi, . . . , Zn] is actually a polynomial. 
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Then notice that for any M G V one still has that M[ip] is holomorphic in z%. 
In particular, M[ifj] is always a polynomial in C(x±, . . . , x n )[zi, ... ,z n ] multiplied by 
the exponential function exp X{Zi. Putting this all together, since we have already 
seen that Q[ip] = q(zi, . . . , z n )ip, we conclude that q(d\, . . . , d n ) G T>q is a constant 
coefficient differential operator. 

Finally, if q G T>q has the property that K o q o G X> then this implies that 
q G Ro(K) and hence that Q G R(K) in the first place. So no operators outside of 
R(K) have the property that they commute with every element of R(K). □ 

2.3. An Explicit Example. It is still not clear from the theorem above whether the 
maximal commutative subrings of T> constructed by Darboux transformation require 
an infinite number of generators over C. Certainly in trivial cases (e.g. K G T> ) the 
resulting ring may require only a finite number of generators. But it would be nice 
to prove that this is always the case or alternatively to observe at least one example 
which does not. By considering a particular example in detail here, we achieve the 
latter. 

2.3.1. A Subring ofC[x,y]. Notation: Let TZ X C C[x,y] (A G C) be the subset 

K\ = {q(x,y) G C[x,y] | q x (z, -) = q y (z, -) = q xy (z, -) = 0}. 

z z z 

In other words, 1Z\ is the set of polynomials q G C[x, y] such that q x , q y and q xy all 
have a factor of xy — A. Note that C C 1Z\ and more importantly that if qi, q 2 G 
are two such polynomials then q\ + q 2 G 1Z\ and q\q 2 G 1Z\. This obviously gives us 
that 

Lemma 2.2. 1Z\ is a proper subring ofC[x,y] containing C as well as every poly- 
nomial of the form p(x, y)(xy — A) 3 for p G C[x, y\. 

It will be shown below that for a particular choice of K G V the maximal com- 
mutative subring R{K) (I'D is isomorphic to 1Z\. Therefore it is interesting to note 
that this ring requires an infinite number of generators over C. 

Lemma 2.3. The ring 1Z\ has the form C[ui(xy — A) 3 , uj 2 {xy — A) 3 , . . . ] where {uji} 
is any basis ofC[x, y] as a vector space. In particular, 1Z\ is not finitely generated. 

Proof. We must show that a polynomial q G C[x, y] is in TZ\ if and only if it is of the 
form g(x, y)(xy — X) 3 +c for some g G C[x, y] and c G C. Clearly, such a q is an element 
of 7Z\. Alternatively, let us suppose that q G 7Z\ and therefore q x = (xy — X)r(x, y). 
Then, since 



also has a factor of xy — A one finds that r has a factor of xy — A and hence q x actually 
has a factor of (xy — A) 2 . (Similarly for q y .) 

Now we have that q x = (xy — \) 2 g(x, y) for some g G C[x, y}. Then integrating by 
parts with respect to x one has 



q xy = xr(x, y) + (xy - X)r y (x, y) 



q(x,y) 
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Continuing to integrate by parts (choosing always to integrate (xy — A) J so that one 
gets higher powers of xy — A and higher derivatives of g) one gets a finite sum (since 
a high enough derivative of g will eventually vanish) of terms each having a factor of 
(xy — A) 3 , plus a constant of integration at the end. 

Now we note that TZ\ cannot be constructed by a finite number of generators over 
C. We may, for instance, suppose that R = C[u\, v 2 , . . . ] for some polynomials z/j and 
w.l.o.g. we may take z/j to have no constant term (z^(0, 0) = 0). But then consider the 
polynomials x l (xy — A) 3 which are elements of These cannot involve any products 
of the generators z/j else they would have a factor of xy — A to a higher degree. Thus, 
they must be a linear combination of some generators v { . On the other hand, since 
the polynomials x l (xy — A) 3 are linearly independent for different i's, it follows that 
you need infinitely many generators to construct 1Z\. □ 

2.3.2. Isomorphism to a ring of Differential Operators. Let us use the notation of 
the preceding subsection to describe a maximal commutative ring of differential op- 
erators. For this example, we will be working in two dimensions only, so n = 2. The 
constant coefficient differential operator which we will factor is p(d± : d 2 ) = (d\d 2 — A) 3 
(A G C) which factors as p = L o K with 

K = xix 2 (did 2 - A) o 



and 



XiX 2 

L = d x 2 d 2 2 + -d ± d 2 2 - x x ~ 2 d 2 2 + -d 1 2 d 2 
x 1 x 2 

I-2AX1X2 -1-Axix 2 „ -2o 2 , -l-Axia^a 
H did 2 H d 2 -x 2 di H — di 

X X X 2 Xi 2 X 2 x x x 2 l 

, 2 1 A 
+A 2 + 



X^X 2 2 X\X 2 



Lemma 2.4. A constant coefficient operator q(d\, . . . ,d n ) G V is an element of 
Ro(K) if and only if the function 



,112+127 



ip(x 1 ,x 2 ,z) := Xix 2 e 2 

is in the kernel of the operator K o q for all values of z G C. 

Proof. One direction is especially simple. HKoq = QoK then 

K o q[ Xl x 2 e xlz+x ^} = Qo x^d^ - X)[e XlZ+X2 ^] = 0. 

Conversely, let us suppose that K o q annihilates this function. This means that 
M := K o q o x\x 2 applied to exp[x\Z\ + £2-22) is zero f° r & H z i z 2 — A = 0. But note 
that M applied to this exponential results in a polynomial in Zi with coefficients in 
C(xi, x 2 ) multiplied by an exponential. This product vanishes on z±z 2 — A = if and 
only if the polynomial has a factor of z±z 2 — A which implies that 

M = Q o xix 2 (9i9 2 - A) 

for some Q G T>. Multiplying this equation on the right by on the right proves 
the lemma. □ 
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Using this lemma and the previous theorem, as well as the bispectrality [fTO] , |TT|| 
of the constant coefficient operators, we demonstrate an isomorphism between R(K) 
and 1Z\. 

Theorem 2.2. The ring R(K), known to be maximal commutative by the preceding 



theorem, is isomorphic to the ringlZ\ (cf. Lemma lOj. 

Proof. Using Lemma |2.4| and Theorem |2.1| , we know that R(K) is isomorphic to the 
ring 

R (K) = {qeV \Ko q[ Xl x 2 e XlZ+x * Xz ] = 0}. 

However, since d Zi := -j^- commutes with differential operators in the variables x,, 
this property is equivalent to saying that 

8 Zl d Z2 [K o q [ e ^+^]] = o \/ Zl z 2 - A = 0. 

This can be written as differential equations for q by applying all of these differential 
operators, clearing the denominator by multiplying by a polynomial in x±,x 2 and 
looking at the coefficients of each monomial in X{. These will be differential expres- 
sions for polynomials in Z\ and z 2 including q which must vanish on Z\Z 2 — A. For 
this to happen, it is necessary and sufficient that q x , q y and q xy all have ziz 2 — A as 
a factor. □ 

3. Resultants of Commuting Differential Operators 

In this section we give a definition of resultants for partial differential operators 
(cf. U) including "spectral parameters" and their significance in the commutative 
case. 

3.f . Definitions. Fix 0<nGN and denote by Vl d the -component vector 

^=(^l>2,---) 

where uf run over all monomial, monic differential operators in the variables 
of degree less than or equal to d. In other words, 

By writing it as a vector, we are supposing that they have an ordering. Such an 
ordering is a choice, not determined canonically, but the particular choice is not 
important to the following. Then, for any differential operator L of order d or less, 
we denote by v^(L) the vector whose i th entry is the coefficient of uf in L. In 
particular, L = Vd{L) ■ Q d . Let L\, . . . , L n+ \ be differential operators in the variables 
xi, . . . ,x n having orders h, . . . , l n+ \ respectively. 

Let N := — n + ^ h and construct the matrix R^ = R^Li, . . . , L n+ i) whose rows 
are v N (uf' k o (L t - m)) for all I < % < n + I and all 1 < j < ( n+I ^~ h ). We call 
any maximal minor determinant of R^ a partial [i-shijted differential resultant. Note 
that each partial /x-shifted differential resultant is a polynomial in the variables fii 
(1 < i < n + 1) with coefficients that may depend on Xj (1 < j < n). We define the \i- 
shifted differential resultant of the operators Li to be the polynomial in the variables 
/ij which is the greatest common divisor of all of these maximal minor determinants. 
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Note: In the case Lj G C[<9i, . . . , d n ], this definition is a special case of the polynomial 



resultant [[UJ and that in one dimension with variable coefficients it reproduces the 
differential resultant of ordinary differential operators used to construct the spectral 
curve |20], [21]. The definition of the differential resultant of the operators Li given 
in |J is, in our terminology a particular partial //-shifted differential resultant of the 
operators Li with all /ij = 0. 

3.2. The Commutative Case. As in the one dimensional case, we will here show 
that the differential resultant provides a polynomial equation satisfied by the op- 
erators Li in the case that they mutually commute. The remainder of the section 
will then be comprised of examples and counter-examples of what we would hope to 
have as a consequence. First, following the approach used in M, we demonstrate the 
following essential lemma: 

Lemma 3.1. Any partial ^-shifted differential resultant of the operators Li, . . . , L n+ i 
can be written as 

n+l 

(2) ^Ao(L 8 -ft) 

i=l 

for some partial differential operators Di with coefficients depending on fij (1 < j < 
n + l) and Xk (1 < k < n). 

Proof. Let j be the integer 1 < j < ( n *^ ) such that = 1 is the differential operator 
of order zero in the vector Or . Construct the matrix M of size 

is the identity matrix except for the fact that the j th column is replaced by the vector 
Q N . Note that detM = 1. 

Let R be a maximal square minor of the matrix R^Li, . . . , L n+ i). Note that the 
elements in the j th column of the matrix R ■ M are all monic monomial differential 
operators composed with the operators Li — fa. Then, expanding down this column 
while taking determinants, one finds exactly something of the form (Q) with the 
coefficients of Di coming from the other minor determinants of R. 

On the other hand, it is an elementary fact of linear algebra that det R = det R ■ M 
and so (|2]) must actually be equal to the order zero operator which is the partial 
/i-shifted differential resultant of the operators Li. □ 

Now suppose that the operators Li (1 < i < n + 1) mutually commute. By defini- 
tion, any partial //-shifted differential resultant of these operators is a polynomial in 
the variables /ij with coefficients possibly depending on the variables Xj. As a conse- 
quence of Lemma [JTT] we then find that the operators satisfy this polynomial. 

Theorem 3.1. Let p(fii, . . . , m n+ i) be any p,- shifted differential resultant of the mu- 
tually commuting operators Li, then 

p(L u . . . ,L n+1 ) = 0. 

Proof. Only commutativity of the /tj's with the L/s is required to rewrite p(/ii, . . . , fJ> n +i) 
in the form (|2J). So, since [Li, Lf\ = we can write p{L\, . . . , L n+ i) by substituting 
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Li for fa in (0). This, however, is clearly zero since every term has a factor of Lj — fa 
for some i. □ 

In the one dimensional case, we can moreover say that the //-shifted differential 
resultant is a polynomial in fa and fa with constant coefficients or a multiple of 
such a polynomial by a function of X\. Here, the results proved thus far leave open 
the possibility that the differential resultant will only produce polynomial equations 
satisfied by the operators with explicit dependence on the variables Xj. We were not 
able to produce any such examples or exclude the possibility. 

3.3. The Zero Possibility. As stated in the introduction, given two commuting 
ordinary differential operators L\ and L 2 the determinant of R^Li, L 2 ) (which hap- 
pens to always be square in the case n = 1) is a non-zero polynomial in fa and fa 
which is satisfied by the operators. Here we will see that the differential resultant 
does not always give such useful information in the higher dimensional case. 
Consider the case n = 2 and 

Li = d\ — d 2 — 1 L 2 = d\ o L 1 L 3 — d 2 o L\. 

Note that these operators satisfy the equation L\ — l\ — L\ — L\ = and so one 
might hope, given Theorem |3.1] , that the differential resultant of these operators is 
lA ~ ~~ A*i — A*i ( or at least is a non-zero multiple of this). 

Lemma 3.2. The differential resultant of the operators Li is the zero polynomial in 
the variables fa (1 < i < 3). 

Proof. One could, of course, merely compute the resultant according to the definition. 
However, there is a more direct and informative way to observe this fact. Since these 
operators are constant coefficient, the problem reduces to a problem of polynomial 
resultants. In particular, the resultant is the same as the resultant of the homogeneous 
polynomials 

Pi(#i, x 2 , x 3 ) —x\ — x\ — (/ii + \)x\ P2(xi, X 2 , Xs) — x\ — X\X% - x\x\ - n 2 x\ 

p 3 (x 1 ,x 2 , x 3 ) = x\x 2 — x\ — x 2 x\ - p 3 x\ 
However, it is well known |16| that this resultant will be zero iff these polynomials 



have a common zero in projective space. Although it is true that no "finite" point 
(with x 3 ^ 0) is a common solution to these polynomials for all values of /ij, there 
are solutions at infinity. In particular, note that the point (1,-1,0) satisfies all three 
polynomials regardless of the values of fa. □ 

It is interesting to note the geometry behind this situation. This problem of having 
a zero resultant never arises in the one dimensional case essentially because only one 
point is being added at infinity and that point is never a solution of the homogeneous 
polynomial. Whereas, in higher dimensions, there is "room" at infinity for many 
solutions. 

Note that the same problem can also occur in a non-constant case (and so not 
simply an example of a polynomial resultant). In particular, the differential resultant 
of any three operators from the ring R{K) described in Subsection 2.3.2| will be zero 



regardless of the values of the variables fa. The mundane explanation of this fact 
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here is merely that the N th powers of di and d 2 never appear in W • li o [L i — /ij) and 
so there are columns of the resultant matrix with all zero entries. 

3.4. Positive Results. A more encouraging example is to consider the operators 

L x = d\ - d\ L 2 = x 2 d 1 + x x d 2 L 3 = L x o L 2 - 1 L 1 7 e C. 

It is a non-obvious fact that [L±, L 2 ] = 0, but given this (which is easily checked) it is 
clear that L 3 also commute and that the three together satisfy a polynomial equation 
p(Li, L 2 , L 3 ) = with 

Then the differential resultant (which can be most easily computed not by finding all 
maximal minor determinants but by the formula D/A where D and A are the minor 
determinants specified in is exactly p 3 (/ii, /i2, ^3). This is very nearly what we 



would want (although there is presently no theory to explain the exponent "3" which 
arises) . 

It is intriguing and surprising that the resultant is independent of the variables 
X\ and x 2 in this case. In the one dimensional case, the resultant of two monic 
differential operators is independent of x if and only if the operators commute. Here, 
the situation involves one operator, L 1; with constant leading coefficients and others 
that are not, which cannot happen in the one dimensional case. 

4. Mathematica Code 

The following code (and assistance using it) is available by writing to the first 
author (kasmanOmath . cof c . edu). It is useful for performing many calculations with 
differential operators, and here we will only give a very brief description of how to 
use it. First set the value of the variable dimen equal to the number of variables you 
will be working with and then input the file containing the code. For instance: 



In[l] := dimen=3 



0ut[l]= 3 



In [2] := «pdo-ak.m 

String dimen already defined. .using present value 3 
READY: Partial Differential Operators in 3 dimensions 

Then, a differential operator can be entered using DX[1] , DX[2] , . . . as the ele- 
mentary differential operators and X[l] , X[2] , . . . as the corresponding variables. 
For example: 



In [3] := L=X[1] DX[2]+X[2] DX[1] 
Out [3]= DX[2] X[l] + DX[1] X[2] 
In [4]:= Q=DX[1] ~2+DX[2] ~2 
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2 2 
Out [4]= DX[1] + DX[2] 

Note that when entering or reading a differential operator in this notation, it 
is always assumed that all differentiation has been taken and functions are on the 
"left" , even if it may not be written this way. (That is, X[l] DX[1]=DX[1] X[l].) 
The non-commutativity is only apparent when multiplying two differential operators. 
Multiplication of differential operators is achieved using the command pdomult: 

In [6] := pdomult [DX[1] ,X[1]] 
Out [6]= DX[1] X[l] + 1 
In [7] : = pdomult [Q , L] 

2 

18:10 > multiplying by DX[1] 

2 

18:10 > multiplying by DX[2] 
Simplifying. . . 



18:10 > Simplifying {3, 0, 0} term. 

18:10 > Simplifying {1, 1, 0} term. 

18:10 > Simplifying {2, 1, 0} term. 

18:10 > Simplifying {1, 2, 0} term. 

18:10 > Simplifying {0, 3, 0} term. 



2 3 3 

Out [7]= 4 DX[1] DX[2] + DX[1] DX[2] X[l] + DX[2] X[l] + DX[1] X [2] + 

2 

> DX[1] DX[2] X[2] 

Similarly, if you want to apply the operator L to the function f you simply say 
pdoapply [L,f ] . 

Finally, given a list of N operators in N—l variables, the matrix whose minor deter- 
minants give the differential resultant can be computed as Diff Result [LI , L2 , . . . , LN] 
(Note that this command automatically subtracts the indeterminate mu[i] from the 
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jth p era tor so that the differential resultant is a polynomial in these variables. The 
command TakeRandomDeterminants [mat ,n] may come in handy as well, since it 
takes the determinant of n randomly chosen maximal minors of the matrix mat. 
Here is the code: 

(* PDO-AK.M Mathematica Code for Differential Operators 
by Alex Kasman, College of Charleston 
kasman@math.cofc.edu *) 

(* Set dimen to be the number of dimensions and then input this 
file. It will be 2 by default *) 

(* Use X[l] ... X [dimen] as variables and DX[1] ... DX [dimen] as the 
corresponding elementary differential operators. *) 

(* Multiply two operators with pdomult [L , Q] *) 

(* Conjugate one operator by another by pdoconj [L,Q] , which yields the 
operator Q L Q~l if this is a differential operator. *) 

(* Diff Result [{LI ,LN}] gives the matrix whose minor determinants 
have a gcd which is the differential resultant of Ll-mu[l] , L2-mu[2] , 
etc. *) 

If [StringMatchQ [ToString [dimen] , "dimen"] , 
dimen=2 ; 

Print ["String dimen undefined. .. setting to ", dimen," by default..."], 
Print ["String dimen already defined. .using present value ", dimen]] 

(* Call the variables X[l] , X [2] , . . .X [dimen] and the differential operators 
DX[1] ,DX[2] , . . .DX[dimen] *) 

zerovect=Table [0 , {i , 1 , dimen}] 

DXpower [vect_] : =Module [{i> , Product [DX [i] " (vect [ [i] ] ) , {i , 1 , Length [vect] }] ] 

DXv [alphaj : =Module [{i> , Table [DX [i] , {i , 1 , Length [alpha] >] . alpha] 

pdocoef [L_,zerovect] :=L/. DX[i_]->0 

pdocoef [L_,vect_] : =Coef f icient [Collect [( (Expand [L] 
/ . DXpower [vect] ->SPACE) / . DX [i_] ->0) , SPACE] , SPACE] 

pdotermmult [coef _ , vect_ ,M_] : =Module [{i , sof ar} , 
verbose ["multiplying by ", DXpower [vect] ] ; 
For [i=l ; sof ar=M, KLength [vect] +1 , i=i+l , 
sof ar=pdopowmult [i , vect [ [i] ] , sof ar] ] ; 
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coef sofar] 

pdotermmult [0 , vect_ ,M_] :=0 

pdopowmult2 [i_ , 1 , M_] : =D [M , X [i] ] +M DX [i] 
pdopowmult [i_ ,0,M_] :=M 

pdopowmult [i_ ,n_ ,M_] : =(pdopowmult2 [i ,n,M] ) 

pdopowmult2 [i_ ,n_ ,M_] : =pdopowmult2 [i ,n-l , pdopowmult 2 [i , 1 ,M] ] 
maketablef or [M_] : =Module [{max, i , j ,k, sof ar} , 

For [i=l , i<=dimen+l , i=i+l ,max [i] =Exponent [Collect [M,DX [i] ] , DX [i] ] ] ; 
For [i=l ; sof ar=Table [j [k] , {k, 1 ,dimen}] , i<dimen+l , i=i+l , 
sofar=Table [sofar , {j [i] , ,max [i] }] ] ; 
Flatten [sofar , dimen-1] ] 

maketablef or [M_ ,K_] : =Module [{max, i , j ,k, sof ar} , 

For [i=l , i<=dimen+l , i=i+l ,max [i] =Exponent [Collect [M , DX [i] ] ,DX [i] ] 
-Exponent [Collect [K,DX[i]] ,DX[i]]] ; 

For [i=l ; sof ar=Table [j [k] , {k, 1 ,dimen}] , i<dimen+l , i=i+l , 
sof ar=Table [sof ar , {j [i] ,0,max[i]}]] ; 
Flatten [sofar , dimen-1] ] 

pdoexp [11_ ,n_] : =Module [{outL, j} , (verbose ["Raising operator to power 
" ,n] ; For [j=0 ; outL=l , j<n, j=j+l , verbose ["pdoexp : power 
" , j ] ; outL=pdomult [11 , outL] ] ; outL) ] 

makeshortlist [L_] : =Module [{i , inlist , outlist} , verbose ["making short 
list"]; outlist={}; inlist=maketablef or [L] ; 
For [i=l , KLength [inlist] +1 , i=i+l , 

If [MatchQ [ToString [pdocoef [L , inlist [ [i] ] ] ] , "0"] .verbose ["not in 
list" , inlist [ [i] ] ] , outlist=Union [outlist , {inlist [ [i] ] }] ; 
verbose ["makeshortlist : includes ", inlist [ [i] ]]]] ; outlist] 

pdomult [L_,M_] :=Module[{i, j ,LM, liszt}, 

liszt=maketablef or [L] ; 

For [i=l ; LM=0 , i<=Length [liszt] , i=i+l , 

LM=LM+pdotermmult [pdocoef [L, liszt [ [i]] ] , liszt [ [i] ] ,M] ] ; 
pdosimp [LM] ] 

pdotermapply [coef _ , vect_ ,M_] : =Module [{i , sofar} , 
verbose ["Applying " ,DXpower [vect] ] ; 
For [i=l ; sof ar=M, KLength [vect] +1 , i=i+l , 
sof ar=pdopowapply [i , vect [ [i] ] , sofar] ] ; 
coef sofar] 



pdotermapply [0 ,vect_ ,M_] :=0 
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pdopowapply [i_ , n_ , M_] : = (D [M , {X [i] , n}] ) 

pdoapply [L_ , M_] : =Module [{i , j , LM , liszt} , 

liszt=maketablef or [L] ; 

For [i=l ; LM=0 , i<=Length [liszt] , i=i+l , 

LM=LM+pdotermapply [pdocoef [L, liszt [[i]]] , liszt [[i]] ,M]] ; 
Simplify [LM] ] 

pdosimp [0] : =0 

pdosimp[L_] : =Module [{liszt , i , sofar} , Print ["Simplifying. . . "] ; 
liszt=maketablef or [L] ; For [i=l ; sof ar=0 , KLength [liszt] +1 , i=i+l , 
If [MatchQ[ToString [pdocoef [L, liszt [[i]]]] ,"0"] , Print ["."] , 
verbose ["Simplifying ", liszt [ [i] ], " term."]; 

sof ar=sof ar+Simplif y [pdocoef [L, liszt [ [i] ] ] ] DXpower [liszt [ [i] ] ] ] ] ; 
sofar] 

pdodisplay [L_] : =Module [{liszt , i} , 

liszt=maketablef or [L] ; 

For [i=Length [liszt] , i>0 , i=i-l , 

pdotermdisplay [pdocoef [L , liszt [ [i] ] ] , liszt [ [i] ] ] ] ] 
pdotermdisplay [0 , vect_] : =(bleh=0) 

pdotermdisplay [coef _,vect_] :=Print["+ (", Simplify [coef] ,") ", DXpower [vect] ] 

pdoconj [P_ ,K_] : =Module [{L, lco ,LK,KP, i , j , liszt , deg, bleh, slv,holdit} , 

liszt=maketablef or [P] ; L=Sum [lco [liszt [ [i] ] ] 

DXpower [liszt [ [i] ] ] , {i , 1 , Length [liszt] }] ; deg=pdodeg [P] ; 

For [i=l , KLength [liszt] +1 , i=i+l , 

If [Sum [liszt [ [i]] [ [j]] ,{j , l.dimen}] >deg-l, bleh [liszt [[i]]]=0; 

L=L/ . lco [liszt [ [i] ] ] ->pdocoef [P , liszt [ [i] ] ] , bleh [liszt [ [i] ] ] =100] ] ; 

Print ["finding KP"] ; KP=pdomult [K,P] ; Print ["finding LK"] ; 

LK=pdomult [L , K] ; brahms=maketablef or [LK] ; 

For [i=Length [liszt] , i>0 , i=i-l , verbose ["solving for 

lco [", liszt [[i]] ,"]"] ; If [bleh [liszt [[i] ]] ==0, verbose [" . . .Already 

set ..."], For [j=Length [brahms] , j>0, j = j — 1 , 

holdit=pdocoef [LK , brahms [ [ j ] ] ] ; 

If [StringMatchQ [ToString [D [holdit , lco [liszt [ [i] ] ] ] ] , "0"] , Print [" . "] , 
slv=Solve [pdocoef [LK, 

brahms [ [j] ] ] ==pdocoef [KP , brahms [ [j] ] ] , lco [liszt [ [i] ]] ] [[1]] ; 
Print ["replacement ->",slv]; j=-5; L=L/.slv; LK=LK/ . slv] ] ] ] ; 
(*Print ["Okay if this is zero ", pdosimp [LK-KP] ]; *) L=pdosimp [L] ] 

pdodeg [L_] : =Module [{transL , i , j , deg , const} , 
transL=L/ . DX [i_] ->const [i] SPACE ; 
deg=Exponent [Collect [transL , SPACE] .SPACE] ; 
deg] 
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pdoTeX [L_] : =Module [{liszt , i} , 

liszt=maketablef or [L] ; 

For [i=Length [liszt] , i>0 , i=i-l , 

pdotermTeX [pdocoef [L , liszt [ [i] ] ] , liszt [ [i] ] ] ] ] 

pdotermTeX [0 , vect_] :=(bleh=0) 

pdotermTeX [l.vectj :=Print["+ " ,TeXForm[DXpouer[vect]]] 
pdotermTeX [coef _,zerovect] :=Print["+ " ,TeXForm[coef ] ] 

pdotermTeX [coef _ , vect_] :=Print["+ " ,TeXForm [Together [coef]] , " ** 
" , TeXForm [DXpower [vect] ] ] 

Print ["READY: Partial Differential Operators in ",dimen," dimensions"] 

genminusone [deg_, liszt_] : =Module [{list2 , i , j ,n}, 
list2={}; 

For [i=l , i<=Length [liszt] , i=i+l , 

n=Sum [liszt [[i]] [[j]] ,{j , 1 .Length [liszt [ [i] ]] }] ; 
list2=Flatten [{list2 , Table [Flatten [{liszt [ [i] ] , j } , 1] , 
{j,0,deg-n}]},l]]; 
list2] 

degcomplete [deg_ , entry_] : =Module [{i} , 
Flatten [{entry , deg-Sum [entry [ [i] ] , 
{i , 1 , Length [entry] }] > , 1] ] 

AllMonomialsDegExactly [deg_] : =Module [{list2 , list3, i , j ,n} , 
If [dimen==l,DX[l] "deg, 
list2=Table [{i} , {i , , deg>] ; 

For [i=l , i<dimen-l , i=i+l , list2=genminusone [deg, list2] ] ; 
Table [degcomplete [deg, list2 [ [i] ] ] , {i , 1 , Length [list2] >] ] ] 

AllMonomialsDegAtMost [deg_] : =Module [{list2,i , j ,n}, 
list2=Table [{i} , {i , , deg}] ; 

For [i=l , Kdimen, i=i+l , list2=genminusone [deg, list2] ] ; 
list2] ; 

multbynextdxupto [deg_ , liszt _] : =Module [{list 2 , list 3, i , j ,k,n,n2,n3,n4, op, rem} , 
list2={}; For [i=l,i<=Length [liszt] ,i=i+l, op=liszt [ [i] ] [ [1] ] ; 
list3=liszt[[i]] [[2]] ; n=Length[list3] ; k=Sum [list3 [ [j] ] ,{j ,l,n}] ; 
rem=op ; list2=Flatten [{list2 , {{rem, Flatten [{list3 , {0}} , 1] }}} , 1] ; 
For [ j = 1 , j <=deg-k , j = j +1 , rem=Simplif y [D [rem , X [n+1] ] +rem*DX [n+1] ] ; 
list2=Flatten [{list2 , {{rem, Flatten [{list3 , { j}} , 1] }}} , 1] ] ] ; list2] 

MultByAllMonomialsUpto [deg_ , L_] : =Module [{i , j , k , list2 , rem} , 
rem={{L,{0}}}; verbose ["MultByAllMonomials j = " , 1 , "/" ,dimen] ; 
For [i=l , i<deg+l , i=i+l , rem=Flatten[ {rem, 
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{{D [rem [ [Length [rem] ] ] [[1]] ,X [1] ] +rem [ [Length [rem] ] ] [[1]] 

DX[1] , rem [ [Length [rem] ] ] [[2]]+l>», 1]] ; For [j=2, j<=dimen, j=j+l, 

verbose ["MultByAllMonomials j=" , j , "/" ,dimen] ; 

rem=multbynextdxupto [deg,rem] ] ; 

Table [Expand [rem [ [k] ] [ [1] ] ] , {k , 1 , Length [rem] >] ] 

pwr [liszt_] : =Product [DX [i] "liszt [ [i] ] , {i , 1 .Length [liszt] }] 

Diff Result [liszt _] : =Module [{i , j ,k,maxdeg,n,n2 ,n3,mm,list2 , list3, list4,rem} , 
n=Length [liszt] ; 

maxdeg=Sum [pdodeg [liszt [ [i] ] ] , {i , 1 ,n}] -dimen; 

mm={>; 

list4={>; 

list2=AHMonomialsDegAtMost [maxdeg] ; 

n2=Length[list2] ; 

For [i=l , i<=n, i=i+l , 

verbose ["Diff result i=" , i , "/" ,n] ; 

list3=MultByAHMonomialsUpto [maxdeg-pdodeg [liszt [ [i] ] ] , liszt [ [i] ] -mu [i] ] ; 

mm=Flatten [{mm, list3} , 1] ] ; 

mm=makeamatrix [mm, list2] ; 

verbose ["The Matrix is:"]; 

If [Length [mm] +Length [mm [ [1] ] ] >30 , 

verbose ["TOO BIG TO PRINT"] , 

verbose [MatrixForm [mm] ] ] ; 

(*TakeAllSquareMinors [mm,list4] *) 

mm] 

makeamatrix [liszt a_ , list2_] : =Module [{liszt , i , j ,k,n, list3, list4,bleh} , 
wid=Length[list2] ; list4=liszta/ .DX [i_] ->0; liszt=liszta-list4; 
list3=Table [DXpower [list2 [ [i] ] ] ->Table [Which [k==i , 1 , True , 0] , 
{k,l,wid>] ,{i, Length [list2] ,2,-1}] ; 

Simplify [ (list4*bleh/ . bleh->Table [Which [k==l , 1 , True , 0] , 
{k,l,wid>] )+(liszt/.list3)]] 

TakeAllSquareMinors [mm_] : =Module [{liszt , 

list2,i, j ,k,n,ll,ww,deter,rem,minlist}, n=0; ll=Length [mm] ; 
ww=Length[mm[[l]]] ; If [ww==ll , {Det [mm] } , If [ww>ll .verbose [" (Note : 
Transposing the matrix for simplicity)"]; 

TakeAllSquareMinors [Transpose [mm] ] , For [k=0,k<=ll-ww,k=k+l , 
verbose ["Doing Young diagrams starting with ",k]; liszt={{k}}; 
For [i=l , i<ww, i=i+l , list2={}; For [j=l,j<=Length [liszt] ,j=j+l, 
list2=Flatten [{list2 , Table [Flatten [{liszt [ [j] ] , 

k>,l] , {k,0, liszt [[j]] [[Length[liszt[[j]]]]]}]>,l]] ; liszt=list2] ; 
verbose ["Computing the corresponding determinants. Listing nonzero 
ones : "] ; For [j=l , j<=Length [liszt] , j=j+l , 

deter=Simplif y [Det [Table [mm [ [11-ww+i-liszt [ [j ] ] [ [i] ] ] ] , {i , 1 , ww}] ] ] ; 
If [StringMatchQ [ToString [deter] , "0"] , , n=n+l; rem[n]=deter; verbose[" 
******************* "] ; verbose [" "] ; verbose ["Nonzero determinant 
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number ",n]; verbose [deter] ]]] ; Table [rem[j] , {j , 1 ,n}] ]] ] 

time :=StringJoin[ToString [Date [] [[4]]] , " : " ,ToString[Date [] [[5]]] ," > "] 

verbose [f_] :=Print [time,f] 
verbose [f _ , g_] : =Print [time , f , g] 
verbose [f _ ,g_ ,h_] : =Print [time ,f , g,h] 
verbose [f _ ,g_ ,h_ , i_] : =Print [time ,f ,g,h, i] 
verbose [f _,g_,h_, i_ , j_] : =Print [time ,f , g,h, i , j] 
verbose [f _,g_,h_, i_ , j_ ,k_] :=Print [time ,f ,g,h, i , j ,k] 

normp [a_ , b_] : =Which [b<a , True , True , False] 

TakeRandomDeterminants [mat_ ,n_] : =Module [{11 , ww , randomyoung , i ,bleh} , 

ll=Length [mat] ; 

ww=Length [mat [ [1] ] ] ; 

verbose [11 , "x" ,ww, " matrix"] ; 

randomyoung :=Sort [Table [Floor [Random []* (11-ww+l)] ,{i,l,ww}] , normp] ; 

For [i=l , i<n+l , i=i+l , 

verbose ["Try number ",i]; 

lst=randomyoung; 

verbose ["Young is ",1st]; 

bleh=Simplify [Det [Table [mat [ [11-ww+i-lst [ [i] ]] ] ,{i,l,ww}]]] ; 
verbose ["Determinant is -> ",bleh]]] 
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